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ABSTRACT 

Planets form in disks around young stars. Interactions with these disks cause them to migrate and 
thus affect their final orbital periods. We suggest that the connection between planets and disks may 
be deeper and involve a symbiotic evolution. By contributing to the outward transport of angular 
momentum, planets promote disk accretion. Here we demónstrate that planets sufñciently massive to 
open gaps could be the primary agents driving disk accretion. Those having masses below the gap 
opening threshold drift inward more rapidly than the disk material and can only play a minor role in 
its accretion. Eccentricity growth during gap formation may involve an even more intímate symbiosis. 
Given a small initial eccentricity, just a fraction of a percent, the orbital eccentricity of a massive planet 
may grow rapidly once a mass in excess of the planet 's mass has been repelled to form a gap around the 
planet's orbit. Then, as the planet's radial cxcursions approach the gap's width, subsequcnt eccentricity 
growth slows so that the planet's orbit continúes to be confincd within the gap. 



1. INTRODUCTION 

Onc of the most important scientific developments of 
the nineties was the discovery of extrasolar planets - plan- 
ets that orbit other stars. Currently, more than a hundred 
are known, along with about a dozen systems containing 
more than one planet. The basic data, masses and orbits, 
have revealed two major surprises: (I) Jovian mass planets 
with short period orbits; (II) Isolated planets with large 
orbital eccentricities. 

Current planet formation theories fail to account for the 
formation of giant planets on short period orbits. Instead, 
there is a general consensus that these planets migratcd in- 
ward from where they were born. Theoretical work done in 
the eighties and nineties established that angular momen- 
tum and energy exchanged between planet and disk at dis- 
crete Lindblad resonances would rcsult in the rapid migra- 
tion of the planet (Goldreich & Tremaine 1980), and that 
almost invariably the migration would be inward (Ward 
1986; Artymowicz 1993b, a). In turn, for sufficiently mas- 
sive planets, this exchange would modify the disk's density 
profile. A familiar aspect of the latter phenomenon is the 
opening of a gap around the orbit of a massive planet 
which locks the orbital evolution of the planet to that of 
the disk, a behavior referred to as type-II migration (Ward 
1997). 

2. ACCRETION 

Disk accretion is driven by the outward transport of 
angular momentum. The mechanism by which this is ac- 
complished in protostellar disks is uncertain. Molecular 
viscosity is far too small to be effective. Magnetorota- 
tional instability is a more plausible candidate, but the 
disk's electrical conductivity may be inadequate to sustain 
it. 

Here we examine the possibility that torques from em- 
bedded planets drive disk accretion. This idea is not new. 
Stimulated by Larson's (Larson 1989) suggestion that spi- 
ral waves might drive disk accretion, Goodman & Rafikov 
(2001) proposed that these waves could be excited by plan- 
ets with masses too small to open gaps. However, they 



ignored the migration of the planets. A general argument 
given below establishes that subcritical planet's would dis- 
appear before significant disk accretion could take place. 

The critical mass required for gap formation in a disk 
without any intrinsic viscosity is less than one earth mass 
(Hourigan & Ward 1984; Ward & Hourigan 1989; Rafikov 
2002). Thus subcritical planets are composed of elements 
heavier than helium which comprise a fraction / of order a 
percent of the disk's mass. Through the torques they exert 
at Lindblad resonances, these planets transfer angular mo- 
mentum outward from material interior to their orbits to 
that external to their orbits. The dominant resonances are 
located about a distance h, the scale height of the disk, in- 
ward and outward from the planet's orbit. Therefore these 
two rates are almost equal. Only the small fractional dif- 
ference, of order h/r, causes the planet to drift inward. 
Thus each planet transfers angular momentum across its 
orbit at a rate that is larger by a factor r/h than that at 
which its angular momentum decays. 

Although the major planet-disk interactions occur at 
distance h from the planet's orbit, density waves carry the 
angular momentum farther away and deposit it at distance 
A where typically h < X < r. So the contribution from 
each planet to the local luminosity of angular momentum 
through the disk is a factor X/h greater than the rate at 
which it loses angular momentum. If subcritical planets 
were the solé source of the disk's angular momentum lumi- 
nosity, the timescale for planet accretion would be related 
to that for disk accretion by a factor F w f{X/h). Good- 
man & Rafikov (2001) estímate X/h ^ 2 for earth size 
planets. Since the total mass of subcritical planets is at 
most a percent of the disk mass (/ < 0.01), they could not 
be responsible for the accretion of more than two percent 
of the disk mass. 

Although subcritical planets are unable to drive disk ac- 
cretion, massive planets might. They reside in gaps and 
acérete at the same rate as the surrounding disk mate- 
rial. If they were the solé source of disk accretion, their 
number density and masses would affect the accretion rate 
but be irrelevant to the final outcome. Suppose there were 
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one giant planct per logarithmic radius interval. Then gap 
widths would be of order the radius and one sided torques 
from individual plancts would have approximate magni- 
tudes given by {Mp/M^fi:r'^{VírY, where Mp and are 
thc masses of the planet and the star, S is the disk's sur- 
face density, Q is the orbital angular velocity, and r is 
the orbital radius. Under these conditions, the angular 
momentum luminosity in the disk would be of order the 
one sided torque. This would yield an accretion timescale 
of order [M^/Mpfü^^ provided Mp < Sr^. Note that 
this timescale is independent of the disk mass and demon- 
strates that Júpiter size planets could cause a disk to ac- 
érete in about a million years. Higher rates would occur 
if the planets were more massive or numerous subject to 
the constraint that the mass in planets does not exceed 
the disk mass. 

3. ECCENTRICITIES AND GAPS 

We adopt standard parameters for use in mimerical eval- 
uations. They are: r = 1 AU, M^/M* = 10^^ for disk to 
star mass ratio where Md is the mass inside r = 1 AU, 
fj, = Mp/Mt, w 10~^ for the planet to star mass ratio, and 
h/r « 0.04 for the disk aspect ratio. Also, in this section, 
we treat the viscosity with the customary approximation 
and choose a = 10~^ which, combined with our disk mass 
and scaleheight, sets a mass accretion rate M « 10~^ so- 
lar masses per year. This cholee is consistent with obser- 
vational determinations of accretion rates onto IMy oíd 
T Tauri stars (Hartmann et al. 1998). Higher accretion 
rates are typical at earlier stages, but any giant planets 
that might have formed then would have been consumed 
by thcir parent stars. As wc are conccrncd with plan- 
ets that survived, it is the later accretion stages that are 
relevant to our investigation . If, as we speculate in the 
previous section, there is no intrinsic viscosity and accre- 
tion is entirely driven by planet torques, then the effective 
a may be even lower than the valué adoptcd abovc. 

Goldreich & Sari (2003) (Hereafter GS) suggest that 
planet-disk interactions might have given rise to the large 
eccentricities of extrasolar planets. The important inter- 
actions are those at Lindblad and corotation resonances. 
For sufhciently small eccentricity, the interactions are lin- 
ear. In this limit, those at ordinary Lindblad resonances, 
which excite eccentricity, are slightly less effective than 
those at corotation resonances, which damp it. However, 
nonlincar saturation of the corotation resonances occurs at 
small eccentricities (Goldreich & Tremaine (1981); Ogilvie 
& Lubow (2003); GS) so a finite amplitude instability lead- 
ing to eccentricity growth is a distinct possibility. Perhaps 
the most serious concern with the scenario by which eccen- 
tricity grows due to planet disk interactions involves thc 
relative importance of eccentricity damping due to apsi- 
dal waves, although it is weaker than previously thought 

GS consider eccentricity cvolution in the context of 
steady-state gaps. Here we argüe that conditions during 
gap formation are more favorable for eccentricity growth. 
This is significant sincc it is the initial stage of eccentricity 
growth that is the most problematic. 

Next we estímate the minimal eccentricity required to 
satúrate the corotation resonances when a gap has just 
formed. At this stage its width, w ~ o/m <C r, is the 



largor of the disk's scale height and the planct's Hill ra- 
dius, each of which are much smaller than the cquilib- 
rium width. Corotation saturation occurs when the den- 
sity gradient is flattened over a scale 5 = (em/x)^/^r. Our 
story involves the comparison of several rates, or inverse 
timescales, namely: 

(i) The rate at which the density gradient is flattened at 
a first order corotation resonance by the corotation torque. 
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(ii) The rate at which viscosity reestablishes the original 
density gradient over a scale 5, 
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(iii) The rate at which a gap of size w is opened by 
principal Lindblad resonances, 

1 dw . / 7- \ 5 
w dt \w 

(iv) The eccentricity growth rate assuming only Lind- 
blad resonances are active, 

1 de , / r \ 4 / M, 



t,L^--^^,^{^yn (3) 



"gap 



t: 



e dt \w 



n. 



(4) 



When corotation resonances are unsaturated, eccentricity 
damps at a rate which is smaller by a factor of about 
twenty (Goldreich & Tremaine 1980). 

For the eccentricity to grow, the corotation resonances 
must satúrate. During gap formation, this requires both 
tsat < tyis and tsat < tgap- Only the former is relevant for 
a steady state gap, and it leads to the criterion derived by 
Goldreich & Tremaine (1981) and Ogilvie & Lubow (2003). 
GS estímate a critical initial eccentricity of about 1% for 
eccentricity growth in equilibrium gaps where torques at 
principal Lindblad resonances balance the large scale vis- 
cous torque. However, as we now show, smaller gaps lead 
to less stringent constraints on the required initial eccen- 
tricity. 

The requirement tsat < tvis is satisfied provided 
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and tsat < tgap provided 
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Thus saturation is most readily achieved at a width 
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provided the initial eccentricity satisfies 
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Saturation of corotation resonances is not a sufñcient 
condition for significant eccentricity growth during gap for- 
mation. In addition, the growth rate of eccentricity must 
exceed that of the gap. This occurs for 
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in other words, once the mass expelled from the gap ex- 
ceeds the planet's mass. Corotation resonances are satu- 
rated at this width providcd the initial eccentricity satisfies 



eo > 
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(10) 
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Significant eccentricity growth requires the initial eccen- 
tricity to satisfy both inequalities (8) and (10). 

Note that eccentricity decays if corotation resonances 
are Icss than 5% saturatcd. However, significant decay 
during gap growth can only occur once w > 2Q[Mp/M¿)r. 
This is avoided provided the initial eccentricity satisfies 
the additional constraint 
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which is relevant only for lighter planets or more mas- 
sive disks (e.g. a planet with mass ten times smaller than 
Jupiter's in a tenth of a solar mass disk). Subsequently, 
wc ignore the possibility of significant eccentricity decay 
during gap formation. 

How does a planet's eccentricity grow as its gap widens 
after the corotation resonances are fully saturated? By di- 
viding the gap growth rate by the eccentricity growth rate 
we arrive at 

/ w - wo \ 

If the disk is ten times more massive than the planet, by 

the time the gap approaches its equilibrium width, compa- 
rable to the radius of the disk, equation (12) predicts that 
the eccentricity would have grown by of order ten e-folds. 
The resulting radial excursions of the planet would cause 
it to pass beyond the gap edges. We expect that before 
this happens, additional damping that wc have not mod- 
eled would come into play and limit eccentricity growth. 
From then on it is likely that the eccentricity would be 
coupled to the gap width such that e ~ w/r. Our con- 
clusions about eccentricity evolution as a function of gap 
width are dcpictcd in figure 1. 

To concludc the discussion of corotation resonances, we 
mention the possibility of resonance overlap. This seems 
to be particularly relevant here as we are interested in 
eccentricity evolution when the gaps are narrow, the frac- 
tional distancc between resonances, w^/r^, is small, and 
the fractional width of the resonances, (/xer/w)^/^, is rel- 
atively large. Resonances overlap if 

e>in^ (13) 
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The evolution of eccentricity in this regime requires fur- 
ther investigation. However, as shown in the figure, for 
our fiducial parameters most of the evolution occurs be- 
low this eccentricity, so resonance overlap probably does 
not play an important role cvcn during gap formation. 

Apsidal wavcs also promote eccentricity damping. GS 
found their influence to be marginal for a gap at its equilib- 
rium width. Since the relative importance of apsidal waves 
varies as the cube of the gap width, they are unlikely to 
play a major role when the gap is well below its equilib- 
rium size. Of course, for a planet to maintain its orbital 



eccentricity, the apsidal waves must also be unimportant 
at the equilibrium gap width. But provided the eccentric- 
ity can grow when the gap is narrow, by the time the gap 
approaches its full width the corotation resonances will 
be fully saturated. Then eccentricity maintenance merely 
requires that eccentricity damping by apsidal waves be 
weaker than eccentricity excitation by first order Lindblad 
resonances, a criterion that is satisfied provided (GS): 
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This requirement is about a factor of twenty less stringent 
than that derived by GS for the apsidal wave torque to 
be less important than the small difference between Lind- 
blad and partially saturated corotation resonance torques. 
In particular, this looser criterion is amply satisfied by 
Júpiter mass planets at r = 1 AU. 

4. PLANET DISK SYMBIOSIS 

We have shown that a symbiotic rclation might cx- 
ist between massive planets and the disks in which they 
form. This relation would explain the similarity between 
timcscalcs estimated to be necded for planet growth and 
those derived from observations of disk accretion, because 
disk accretion would only commence after massive plan- 
ets form. In this scenario, most planets commit suicide 
by promoting the accretion of the disk to which they are 
lockcd. That somc planets survive implies that not all 
of the disk material is accreted. Some other mechanism, 
such as evaporation, might remove the final remnants of 
the disks. 

The interplay between eccentricity excitation and gap 
formation is a more subtle aspect of this symbiosis. Once 
a planet attains a suíficient mass, perhaps a few earth 
masses, it rapidly accretes an envelope of hydrogen and 
hclium and begins to open a gap in its disk. Provided the 
planet's orbit is endowed with a small initial eccentricity, 
the corotation resonances satúrate while the gap width is 
not much larger than the disk scale height h. However, de- 
spite their saturation, significant eccentricity growth does 
not occur until the mass cleared from the gap bccomes 
comparable to the planet's mass. Subsequently, the rate 
of fractional eccentricity increase exceeds that of the gap 
width. From then on it is plausible that the eccentric- 
ity maintains a valué of order the fractional gap width, 
e « w/r. A schcmatic dcscription of eccentricity and gap 
growth is displayed in the figure. 

Eccentricity growth during gap opening alleviates two of 
the major concerns raised by Goldreich & Sari (2003) in 
their discussion of possible eccentricity growth for planets 
in stcady-statc gaps. It reduces both the relative impor- 
tance of the apsidal torque and the required valué of the 
initial eccentricity. The case for eccentricity growth due 
to planet disk interactions becomes more promising with 
this new analysis. 

This research was supportcd in part by an NSF grant 
AST-0098301 and NASA grant NAG5-12037. 
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FiG. 1. — Evolution of a planet's eccentricity versus fractional gap width for a variety of initial valúes of eccentricity. Fast eccentricity growth 
occurs to the right of the vertical Une, where the mass expelled from the gap exceeds the planet's mass, and above the two solid diagonal 
Unes (equations (5) and (6)), where the corotation resonances are fuUy saturated. In that región of the (e,w) plañe, shown as unshaded, 
eccentricity grows exponentially until it approaches the fractional gap width shown by the dotted Une marked e = w/r. The green circle marks 
the criterion for eccentricity growth derived by GS for steady-state gaps. However, since signiñcant eccentricity growth can occur before the 
gap reaches its máximum width, the actual criterion is less stringent, and is indicated by the green square. The criterion for resonance overlap 
is shown as dotted line marked "overlap" . 
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